A hazard of open quantum dynamics: Markov approximations encounter map domains 
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A Markov approximation in open quantum dynamics can give unphysical results when a map acts on a state 
that is not in its domain. This is examined here in a simple example, an open quantum dynamics for one qubit 
in a system of two interacting qubits, for which the map domains have been described quite completely. A 
time interval is split into two parts and the map from the exact dynamics for the entire interval is replaced by 
the conjunction of that same map for both parts. If there is any correlation between the two qubits, unphysical 
results can appear as soon as the map conjunction is used, even for infinitesimal times. If the map is repeated 
an unlimited number of times, every state is at risk of being taken outside the bounds of physical meaning. 
Treatment by slippage of initial conditions is discussed. 
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I. INTRODUCTION 

A Markov approximation assumes that the way a state 
changes in time does not depend on the time when the change 
begins. In open quantum dynamics, this assumption does not 
hold when there are coiTelations between the subsystem be- 
ing considered and the rest of the larger system, because the 
changes of states in time are described by maps that depend 
on the correlations, and the correlations change in time. The 
maps generally are not completely positive and act on limited 
domains 1 1,^. Maps for changes that begin at different times 
can have different domains. When a Markov approximation 
requires the same map to act at different times, it may require 
the map to act on a state that is not in its domain. This can 
give unphysical results. A density matrix for a state can be 
mapped to a matrix that is not a density matrix and can not 
represent a physical state. 

We look at this here in a simple example, an open quantum 
dynamics for one qubit in a system of two interacting qubits, 
for which the map domains have been described quite com- 
pletely 111, 0]. We consider only one particular hazard that 
Markov approximations encounter in open quantum dynam- 
ics. We examine only the Markov property that the same map 
is used at different times. We use maps from the exact dy- 
namics. They are not approximations. We just use them out 
of place. We split a time interval into two parts and replace the 
map for the entire interval with the conjunction of that same 
map for both parts. If there is any correlation at all between 
the two qubits, unphysical results can appear as soon as the 
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map conjunction is used, even for infinitesimal times. If the 
map is repeated an unlimited number of times, every state is 
at risk of being taken outside the bounds of physical meaning. 
In this simple example, everything that happens can be seen 
clearly. 

A Markov approximation is not expected to work very well 
for this example because the reservoir is no bigger than the 
system and changes as fast as the system does. There can be 
large congelations between the system and the reservoir and 
they can change just as fast. The exact dynamics of the system 
does not depend on time as a semigroup. Our interest is not in 
a failure of accuracy of a Markov approximation. It is in the 
reason for unphysical results. We show that unphysical results 
can occur simply because a map is made to act on a state that 
is not in its limited domain. Our focus is on map domains. We 
use this two-qubit example because it is the only example we 
have where the domains are known. 

Markov approximations in open quantum dynamics en- 
counter many hazards. Knowing the nature of this one partic- 
ular hazard, seen clearly in this example, may make it easier 
to navigate around it in other settings. Separating it from other 
hazards may make it possible to view them more clearly. 

To be able to see what happens, we use maps of mean values 
to describe both the open dynamics of the single qubit and the 
relevant parts of the full dynamics of the two qubits. This 
keeps us within sight of established navigation marks. The 
map domains are described in terms of mean values. They 
depend on coiTelation mean values that make the maps change 
with time. These marks are not visible when you work with a 
master equation or use an operator sum form to describe the 
maps. 

The importance of correlations has been long appreciated 
iHQ]- The role of map domains has not been described. One 
measure used to deal with hazards of unphysical results is to 
adjust or "slip" the initial conditions IH 0, [M 1^ • This is 
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discussed, in the light of the example, in Section IV. First the 
example is set out, the exact dynamics in Section II, and the 
map conjunctions in Section III. 

II. EXACT DYNAMICS 

We consider two qubits, one described by Pauli matrices 
Si, S2, S3 and the other by Pauli matrices Si, S2, S3, and 
consider the dynamics generated by the Hamiltonian 

H = ^EsEi. (2.1) 

We focus on the open dynamics of the S qubit described by 
maps of the mean values (Si), (S2), (S3) at time to 

(51) (t) = (Si)cosi- (S2Si)sint 

(52) (i) = (S2)cosi+ (SiSi)sint 

(53) (i) = (S3) (2.2) 

at times t. We write (S) for the vector with components (Si), 
(S2), (S3), or for the set of those three components, and write 
(S)(i) for (Si)(<), (S2)(t), (S3)(t). For each t and each 
(S2S1) and (SiSi), there is a map from states at time 0, de- 
scribed by mean values (S), to the states at time t described 
by the mean values (S)(i). Different (S2S1) and (SiSi) give 
different maps. 

When there are correlations between the two qubits, these 
maps generally are not completely positive and act in limited 
domains. Each map is made to be used for a particular set 
of states described by a particular set of (S), which we call 
the compatibility domain. It is the set of (S) that are com- 
patible with the (S2S1) and (SiSi) in describing a possible 
initial state for the two qubits. In a larger domain, which we 
call the positivity domain, the map takes every positive matrix 
to a positive matrix. The positivity domain is the set of (S) 
for which |(S)(t)| < 1. We have described these domains 
quite completely for the maps we consider here [Jj, 2j\- For 
these maps, the compatibility domain is the intersection of all 
the positivity domains for different t for the same values of 
(S2S1) and (SiSi); this is not generally true |[lll2l- 

Here, with no loss of generality, we let (S2S1) be zero and 
(SiSi) positive. Then the maps depend on just the one pa- 
rameter (Si Si). It changes in time to 

(SiSi)(i) = (SiSi) cost - (S2) sint. (2.3) 

III. MAP CONJUNCTIONS 

Let (S)(<|s) be the replacement for (S)(i + s) obtained 
by going from time to time t with the map established at 
time with (Si Si) and then from time t to time t + s with 
the same map, with (Si Si) unchanged. We call this a map 
conjunction. We use it to examine one property of Markov 
approximations that can give unphysical results: application 
of the same map at different times. This replacement is not a 



Markov approximation itself. The semigroup property is used 
only at the point where the maps are joined. The maps are 
from the exact dynamics. They are not approximations. They 
are just being used out of place. 

We look at this map conjunction for initial states where 
(Si) and (S3) are zero. Then (Si)(t), (S3)(t) and (Si)(t|s), 
(S3)(t|s) are all zero, and 

{^2){t\s) ^ (S2)(t) COS s + (SiSi) sins 
= (S2) cost cos s 

+ (EiSi)(sint coss + sins). (3.1) 

When (El) and (S3) are zero, (S) is in the compatibility 
domain for (SiSi) if IBll 

(S2)2 + (SiSi)2 < 1. (3.2) 

To push to the limit, we look at a point on the edge of the 
compatibility domain and let 

(SiSi) = sing, (S2)=cos(7, (3.3) 

both positive. In particular, we look at what happens when t 
is q. We have 

(S2)(g) = 1 

(SiSi)(g) = 
(S2)(g + s) = coss 

(E2)(g|s) = coss + sing sins. (3.4) 

When s is 0, both {T.2){q + s) and {T,2){q\s) are 1. Then, 
as s increases, (S2)(<7 + s) goes down, but (S2)(g|s) goes 
up; when s is 0, the slope d{J^2){q\s) /ds is sing, which we 
are assuming is positive. We have the physically impossible 
result that {T,2){q\s) becomes larger than 1; it can not be a 
mean value (S2) for a physical state. This happens for any 
(SiSi) that is not zero. It happens for infinitesimal s. If there 
is any correlation at all between the two qubits, unphysical 
results appear as soon as the map conjunction is used. 

The unphysical results occur when the map acts on a state 
that is not in its domain. From Eqs.( |3.2| ), we see that when t is 
q and (S2) (t) is 1, the state is not in the compatibility domain 
of the map. It is in the compatibility domain for the map that 
could be established at time q to follow the exact dynamics, 
because (SiSi)(t) is when t is q. 

Now we look at what can happen with more repeti- 
tions of the map. Let (S)(i|si|s2) be the replacement for 
(S)(<; + si + S2) obtained by going through three steps in 
time from Oto<to< + si toi + Si + S2 using the same 
map established at time 0, with the same (SiSi), for each 
step. In general, let (S)(t|si| ... |s„) be the replacement for 
(S)(t + si + ... + s„) obtained by going through n + 1 steps 
in time from to < to i + si ... to t + si + ... + s„ using that 
same map for each step. Now we want to look at states deep 
inside the compatibility domain, not on the edge, so we do not 
assume Eqs. ( 13.31 ), and we do not assume that t is q. We do 
still assume that (Si) and (S3) are zero. From Eqs. ( 12.21 ). we 
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can see that as t varies, the magnitude of (S2)(t) is maximum 
when 

[(I]2)W]' = (S2)2 + (EiSi)2. (3.5) 

From this and the first Hne of Eq. (13. Il l, we can see that as t 
and s vary, the magnitude of {J^2){t\s) is maximum when 

= (S2)2 + 2(I]iSi}2. (3.6) 

We can see similarly that as t, si and §2 vary, the magnitude 
of {Y,2){t\si\s2) is largest when 

m){t\s,\s2)r - [(S2)(iki)]' + (SiSi)2 

- (S2)2+3(EiSi)2, (3.7) 

and as t, si, ... and s„ vary, the magnitude of 
{T,2){t\si I ... \sn) is largest when 

[(E2)(iki|...|s„)]2 = [(E2)(iki|...|s„_i)]2 + (EiSi)2 
= (S2)' + (n+l)(I]iSi)2. (3.8) 

As n increases, the magnitude of (E2)(i|si| ••• \sn) can grow 
with each step and eventually become larger than 1. This can 
happen no matter how small (E2} and (EiSi) maybe, as long 
as (EiSi) is not zero. If there is any correlation at all between 
the two qubits, every state of the S qubit is at risk of being 
taken outside the bounds of physical meaning. 

Again, unphysical results appear when the map acts on a 
state that is not in its domain. From Eqs. ( I3.2l i and ( I3.8l l, 
we can see that when the magnitude of {J^2){t\si\ . . . |s„) is 
larger than 1, the predecessor (I]2)(^|si| ■ • ■ |sn-i) is not a 
mean value (S2) for a state in the compatibility domain of 
the map. Repetitions of the map can move a state out of the 
compatibility domain from deep inside it. 

IV. SLIPPING INITIAL CONDITIONS 

One measure used to deal with hazards of unphysical re- 
sults is to adjust or "slip" the initial conditions, or restrict 



what they can be fi |^ 0, S, 1^ ■ The maps are kept the same. 
In our example, we see that unphysical results can indeed be 
eliminated by changing or restricting the initial value (S2). 
From Eqs. ( 13.21) . (13.5b and ( 13.61 ). we see that the magnitude of 
(E2) remains smaller than 1 for all t and s if 

(S2)2 + 2(EiSi)2 < 1, (4.1) 

which means that the state represented by is in the 

compatibility domain. This is more restrictive than the re- 
quirement (I3.2I 1 that the state represented by (S) is in the com- 
patibility domain, but it may be an acceptable restriction when 
the correlation between the two qubits described by (EiSi) is 
small. From Eq.( l3.8l l. we see that 

(S2)2 + (n+l)(SiSi)2 < 1 (4.2) 
would be needed to prevent unphysical results when n repeti- 
tions of the map are arranged to occur at times when they do 
maximum damage. Such an arrangement is not likely to occur 
in a well constructed model, but knowing how it could occur 
allows care to be taken to see that it does not. 

In cases where unphysical results have been eliminated by 
restricting initial conditions, the maps of the remaining initial 
states have been observed to produce entanglement with an- 
other system that is separate from the dynamics i^l^. This 
is not necessarily a symptom of pathology. Such entangle- 
ment can be produced by the maps that describe the correct 
exact open quantum dynamics of a subsystem when there are 
correlations between that subsystem and the rest of the larger 
system involved in the dynamics lilOll . 
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